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5_l ' In this talk we discuss at hand of two examples the crucial role played by self consistency 

f> || in hadron physics. The first example concerns the quark-mass dependence of the baryon 

. octet and decuplet masses. It is shown that within a self consistent one-loop approach based 

on the chiral Lagrangian the 'mysterious' quark-mass dependence of the S mass predicted by 
' the MILC collaboration may be recovered in terms of a discontinuous chiral extrapolation. 

This is a consequence of self consistency imposed on the partial summation, i.e. the masses 
used in the loop functions are identical to those obtained from the baryon self energies. In 
the second example we discuss recent studies on the properties of D mesons in cold nuclear 
matter as they are predicted by coupled-channel dynamics. Here a self consistent many-body 
approach reveals the close interlink of the properties of D meson and open-charm baryon 
Qjj ■ resonances in nuclear matter. The relevance of exotic baryon resonances for the spectral 

^ ' distortion of the in nuclear matter is pointed out. 

§1. Introduction 

In this talk we discuss two seemingly unrelated topics. First the quark-mass 
dependence of baryon masses and second the properties of D mesons in cold nuclear 
matter. The two topics exemplify the crucial importance of self consistency in hadron 
physics. As it will become clear during the talk both topic have exotic aspects 
■ justifying its presentation in a workshop about 'exotics', 

t*"^ | The present-day interpretation of QCD lattice simulations requires a profound 

understanding of the dependence of observable quantities on the light quark masses. 
A powerful tool to derive such dependencies is the chiral Lagrangian, an effective 
field theory based on symmetry properties of QCD. The application of strict chiral 
perturbation theory to the SU(3) flavor sector of QCD is plagued by poor convergence 
properties for processes involving baryons. 1 - 1-3 ) Thus it is important to establish 
partial summation schemes that enjoy improved convergence properties and that are 
better suited for chiral extrapolations of lattice simulations. We review the results 
of 4 ). 5 ) where recent lattice QCD simulation of the MILC collaboration, 6 )' 7 ) that 
use dynamical u-,d- and s-quarks in the staggered approximation are interpreted. 
Adjusting the values of the Q 2 counter terms to the physical masses we suggest 
the possibility of a discontinuous dependence of the baryon masses on the pion 
mass. This is a consequence of self consistency imposed on the partial summation 
approach, i.e. the masses used in the loop function are identical to those obtained 
from the baryon self energy. The latter is a crucial requirement since the loop 
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functions depend sensitively on the precise values of the baryon masses. Our results 
may explain the mysterious quark-mass dependence for the E mass observed by the 
MILC collaboration. 6 )' 7 ) 

Except for the recent works 8 ) -10 ) a mean field ansatz was assumed to predict 
the properties of D mesons in nuclear matter. 8 )' n )~ 13 ) This development resembles 
to a large extent the first attempts to predict the properties of kaons and antikaons 
in nuclear matter (see e.g. 14 ) -16 )), which also assumed a mean- field type behavior 
for the mass shifts in nuclear matter at first. However, by now it is well established 
that for the antikaon such an ansatz is not valid due to complicated many-body 
dynamics induced by the presence of resonance-hole states. 17 )' 18 ) A self consistent 
approach is required that treats the spectral distortions of the hyperon resonances 
and the antikaon in a consistent manner. As illustrated in 9 ) it is crucial to study also 
the properties of D mesons in a more microscopic manner. The key ingredient of a 
realistic description are the D-meson nucleon scattering amplitudes, that determine 
the D meson self energy at least for dilute nuclear matter. 16 )' 19 ) In contrast to the 
D_-nucleon, the D + -nucleon system may couple to open-charm baryon resonances. 
This poses a particular challenge since the spectrum of the open-charm baryon res- 
onances is studied experimentally and theoretically poorly so far. As suggested first 
in 20 ) the well established /l c (2594) s-wave resonance resonance may be dynamically 
generated by coupled-channel interactions. 8 ) A detailed study 21 ) that took into ac- 
count charm-exchange reactions systematically revealed that it couples strongly to 
the .D_-nucleon channel. Thus the /l c (2594) should play an important role in the 
description of the nuclear D_ dynamics. Besides the narrow yl c (2594), which couples 
dominantly to the DN, D S A channels, the work 21 ) recovers a broad s-wave state, so 
far unobserved, that is interpreted as a chiral excitation of the open-charm sextet 
ground states. 20 ) The latter couples strongly to the ttE c channels but very weakly 
to the DN channel. It is the analogue of the yl(1405), which is a chiral excitation of 
the baryon octet ground states. 22 ) -25 ) A further striking prediction of the work 21 ) is 
a narrow isospin one resonance of mass 2.62 GeV which couples dominantly to the 
DN,D S E channels. This state is so far unobserved, but, so it existed, would affect 
the properties of D + mesons in nuclear matter significantly. 9 ) 

§2. Chiral extrapolation of baryon masses 

We collect the terms of the chiral Lagrangian that determine the leading orders 
of baryon octet and decuplet self energies. 26 ) Up to chiral order Q 2 the baryon 
propagators follow from 

C = tr (B [i$ - M [8] ] B) 

- tr (X • ( [i - M m \gT~i (7 M ^ + + ^ [i $>+ M [w] ] Y) 4,) 

- 2 do tr (X • A^ tr ( Xo ) - 2 d D tr ((4, • A") X o) 

+ 2 6 tr(5 5) tr(xo) + 2 6 F tr(s [ Xo , B]^j + 2 b D tr (b { X o, £}) , 
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I m 2 J . (2-1) 

2m\-m\ 

where we use the notations of 4 )' 5 ) for the baryon octet and decuplet fields B and 
A respectively. The evaluation of the baryon self energies to order Q 3 probes the 
meson-baryon vertices 

C = ^ tr (B 757^ [d^, B]) + ^ Tr (b 75 7 M {<9 m #, B}) 
- tr • {d v $) [sT -\z-f-f]B + h.c. 

-^tr([^- 7 57,^](^)), (2-2) 

where we apply the notations of. 27 ) We use / = 92.4 MeV in this work. The values 
of the coupling constants F, D, C and H may be correlated by a large- iV c operator 
analysis. 28 )~ 30 ) At leading order the coupling constants can be expressed in terms of 
F and D only. We employ the values for F and D as suggested in. 27 )' 31 ) All together 
we use 

F = 0.45, L> = 0.80, H = 9F-3D, C = 2D. (2-3) 

We take the parameter Z = 0.72 from a detailed coupled-channel study of meson- 
baryon scattering that was based on the chiral SU(3) Lagrangian. 27 ) The latter 
parameter is an observable quantity within the chiral SU(3) approach: it contributes 
at order Q 2 to the meson-baryon scattering amplitudes and cannot be absorbed into 
the available Q 2 counter terms. 27 ) 

It is straight forward to evaluate the one-loop contributions to the baryon self 
energy as implied by (|2-2p . Explicit and complete expressions are provided for the 
first time in. 4 )' 5 ) Since we do not advocate a strict chiral expansion, rather a partial 
summation, the renormalization needs to be discussed briefly As pointed out in, 4 ) 
within dimensional regularization there persists an ambiguity on how to implement 
the chiral counting rules. This leads to the presence of an infrared renormalization 
scale that can be used to optimize the speed of convergence. Performing a strict 
chiral expansion the physical parameters are independent on the infrared scale, as 
they should be. However, the size of the Q 2 counter terms depend linearly on this 
scale. In turn the apparent convergence properties reflect the choice of that scale. 
Only for reasonable values of the infrared scale, fim, the counter terms have a natural 
size. 

Once a partial summation, as implied by a self consistent evaluation of the loop 
function, is performed an infinite number of counter terms would be needed to arrive 
at results that are manifestly independent on the renormalization scales. Dialing a 
value of the infrared scale amounts to a particular choice thereof. Any residual 
dependence of the baryon masses on the infrared scale may be used to estimate the 
uncertainty of the given truncation. 
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/x/ij = 350 MeV 


fJ-IR z 


= 450 MeV 


Hm = 550 MeV 


b [GeV- 1 ] 


-0.89 






-0.63 




-0.38 




b D [GeV- 1 ] 


+0.29 






+0.19 




+0.10 




b F [GeV- 1 ] 


-0.34 






-0.25 




-0.15 




do [GeV- 1 ] 


-0.22 






-0.15 




-0.08 




d D [GeV" 1 ] 


-0.35 






-0.30 




-0.24 




M N [MeV] 


750 + 310 - 


121 


813 + 


232 - 


105 


875 + 153 - 


89 




= 939 






= 939 




= 939 




M A [MeV] 


750 + 536 - 


150 


813 + 


398 - 


79 


875 + 260 - 


9 




= 1136 






= 1131 




= 1126 




M E [MeV] 


750 + 882 - 


425 


813 + 


630 - 


239 


875 + 379 - 


54 




= 1207 






= 1203 




= 1200 




M s [MeV] 


750 + 934 - 


361 


813 + 


680 - 


171 


875 + 427 + 


19 




= 1323 






= 1321 




= 1320 




M A [MeV] 


1082 + 241 - 


91 


1108 + 


164 - 


39 


1133+ 87 + 


12 




= 1232 






= 1232 




= 1232 




M s [MeV] 


1082 + 347- 


49 


1108 + 


253 + 


17 


1133+ 160 + 


83 




= 1380 






= 1378 




= 1376 




Ms [MeV] 


1082 + 453 - 


4 


1108 + 


343 + 


78 


1133+ 233 + 


160 




= 1530 






= 1528 




= 1526 




M fi [MeV] 


1082 + 558 + 


34 


1108 + 


432 + 


134 


1133+ 305 + 


235 




= 1674 






= 1674 




= 1674 




cr^jv [MeV] 


52.4 


53.9 


55.7 


[MeV] 


384.0 


380.1 


386.3 


a K - n [MeV] 


359.4 


354.8 


361.5 



Table I. The parameters are fitted to reproduce the baryon masses at physical pion masses as well 
as the SU(3) limit values M [8] ~ 1575 MeV and M [10] ~ 1710 MeV at ~ 690 MeV. We use 
fJ-Uv = 800 MeV. The masses are decomposed into their chiral moments. 

For given values of the infrared and ultraviolet renormalization scales the pa- 
rameters 6d,6f and do are fitted to the mass differences of the octet states and 
decuplet states. The absolute mass scale of the octet and decuplet states can be 
reproduced by appropriate values of the bare baryon masses. This procedure leaves 
undetermined the two parameters bo and do. Good representations of the physical 
baryon masses can be obtained for a wide range of the latter. The parameter bo may 
be used to dial a given pion-nucleon sigma term at physical pion masses. Similarly 
the unknown parameter do may be determined to reproduce a given pion-delta sigma 
term. The latter two parameters are adjusted as to reproduce the baryon octet and 
decuplet masses in the SU(3) limit at m n ~ 690 MeV. The MILC simulations suggest 
the values Mr 8 i ~ 1575 MeV and M[ 10 ] — 1710 MeV. This procedure is biased to the 
extent that it assumes that the self consistent one-loop results will be applicable at 
such high quarks masses. On the other hand as long as there are no continuum limit 
results of the MILC collaboration available this is an economical way to minimize 
the influence of lattice size effects. The latter are expected to be smaller at large 
quark masses. The procedure may be justified in retrospect if it turns out that the 
extrapolation recovers the behavior predicted by the lattice simulation. 

The parameters used are collected in Tab. U together with the implied masses of 
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Fig. 1. The pion mass dependence of baryon octet and decuplet masses predicted by the chiral loop 
expansion taking the parameters of Tab. U The lines represent the masses for the infrared scale 
put at = 450 MeV. The solid squares are the simulation points of the MILC collaboration. 



the baryon octet and decuplet states. A fair representation of the physical baryon 
masses is obtained. As discussed in detail in 4 ) the parameters &o,d,f and <io,D show 
a strong dependence on the infrared scale hir. In contrast, the physical baryon 
masses suffer from a weak dependence only. For natural values of the infrared scale 
the chiral expansion appears well converging as indicated by the decomposition of 
the baryon masses into their moments. Taking the residual scale dependence of the 
pion-nucleon sigma term as a naive error estimate we obtain a^N = 54 ± 2 MeV. 

We turn to the pion-mass dependence of the baryon octet and decuplet masses. 
It should be emphasized that the baryon masses are a solution of a set of coupled and 
non-linear equations in the present scheme. This is a consequence of self consistency 
imposed on the partial summation approach. As a consequence of the non-linearity 
for a given parameter set there is neither a guarantee for a unique solution to exist, 
nor that solutions found are continuous in the quark masses. Indeed, as illustrated 
by Fig. 1 the pion-mass dependence predicted by the chiral loop expansion is quite 
non-trivial exhibiting various discontinuities. It is striking to see that we reproduce 
the 'mysterious' pion-mass dependence of the E mass, i.e. the quite flat behavior 
which does not seem to smoothly approach the physical mass. Given the present 
uncertainties from finite lattice spacing, the staggered approximation and the theo- 
retical uncertainties implied by higher order contributions, we would argue that we 
arrive at a fair representation of the lattice simulation points for all baryons (see 5 )) 
with some reservation concerning the nucleon. Incorporating the many Q 4 counter 
terms offered by the chiral Lagrangian it is reasonable to expect that the latter will 
further improve the picture. However, as long as there is no detailed analysis avail- 
able that performs the continuum limit there is not much point considering the Q 4 
counter terms. 
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§3. D mesons in nuclear matter 

After a demonstration of the crucial importance of self consistency for the chiral 
extrapolation of the baryon octet and decuplet masses we turn to the propagation of 
D mesons in a cold nuclear environment. We discuss the properties of the and 
Df mesons as derived from a self consistent many body approach based on coupled- 
channel dynamics. 9 - 1 The dominant interaction was modelled by the exchange of light 
vector mesons in the t-channel. 21 ' All relevant coupling constants were obtained from 
chiral and large-A'c properties of QCD. Less relevant three-point vertices related 
to the t-channel forces induced by the exchange of charmed vector mesons were 
estimated by a flavor SU(4) ansatz. 21 ) The resulting s-wave D-N — > D^N scattering 
amplitudes are dominated by the dynamically generated 7l c (2594) and I7 C (2620) 
resonances. The s-wave DjrN — > D^N scattering amplitudes are characterized by so 
far unobserved exotic resonances at 2.89 GeV and 2.78 GeV. Only the D_N — > D-N 
scattering process is not influenced by the presence of a resonance. In this case the 
amplitude is characterized to a large extent by the scattering length 

45rJ5 ~ -0.16 fm, a^I^ -0.26 fm. (3-1) 

The many-body computations are based on the self consistent and covariant 
formalism established in. 17 )' 18 - 1 It is important to perform such computations in a 
self consistent manner since the feedback of an altered meson spectral function on 
the resonance structure is typically an important effect. 17 ^' 18 )' 32 ) Numerical results 
for the spectral distributions of the D mesons as derived in 9 ) are shown in Fig. [2j 
Given the scattering lengths (|3-1|) the mass shift of a D_ meson in nuclear matter 




Fig. 2. Spectral distributions of the D and Df 1 mesons. Results are shown for the meson momenta 
0, 200 and 400 MeV and nuclear densities 0.17 fm -3 and 0.34 fm~ 3 . The self consistent many- 
body approach of 18 ' was applied. 
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Fig. 3. Spectral distributions of the D ± and D^ 1 mesons. Imaginary part of the isospin zero (l.h.p.) 
and isospin one (r.h.p) D+-nucleon scattering amplitude at saturation density as compared to 
the free-space case (dotted lines). The amplitudes are shown for two values of the resonance 
three-momentum w — GeV and w = 0.4 GeV. 

is fully determined by the low-density theorem. 16 )> 19 ) The low-density mass shift of 
17 MeV is quite close to the self consistent result shown in Fig. [2j Self-consistency 
leads to a mass shift of 18 MeV at saturation density and 38 MeV at twice saturation 
density. The spectral function of the Z)+ meson has a two-mode structure, which is 
a consequence of important resonance-hole contributions. At saturation density the 
main mode is pushed up by about 32 MeV as compared to the free-space meson. 
Our results for the D + meson differ from the previous study 8 ) significantly. This 
is a consequence of the quite different interaction used in the two computations. 
The work 8 ) overestimates the charm-exchange channels. In 21 ) those channels are 
suppressed by a kinematical factor m 2 p /m 2 D ~ 0.2. In particular the work 8 ) did not 
predict the isospin one resonance I7 C (2620). The latter dominates the resonance-hole 
component in the spectral distribution. 

We turn to the mesons with non-zero strangeness. The study of 21 ) predicts 
the existence of a coupled-channel molecule in the (D~N, DA, DS) system. The 
latter carries exotic quantum numbers that can not be arranged by three quarks 
only. Exotic s-wave states with C=-l were discussed first by Gignoux, Silvestre-Brac 
and Richard 33 ) and later by Lipkin. 34 ) The binding of about 190 MeV for such a 
the state as predicted in 21 ) awaits experimental confirmation. If confirmed the s- 
wave D~N — > D~N scattering amplitude must show a prominent pole structure at 
subthreshold energies. The latter leads to a well separated two-mode structure of 
the Dj spectral distribution. The main mode is pushed up by less than 10 MeV at 
nuclear saturation density. 

The in-medium spectral distribution of the is most striking. The additional 
mode expected from the possible existence of the exotic state at 2.89 GeV merges 
with the main Df mode into one broad structure, in particular at intermediate meson 
momenta 400 MeV. The results are quite analogous to the spectral distribution of 
the K_ where the yl(1405) nucleon-hole state gives rise to a broad distribution. 17 )' 18 ) 
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Clearly this result is an immediate consequence of the small binding energy of the 
exotic state. 

In Fig. [3] the spectral functions of the /l c (2594) and Z" c (2620) resonances are 
shown at saturation density as compared to their free-space distributions. We observe 
small attractive shifts in their mass distributions and significant broadening, at least 
once the resonances move relative to the matter bulk. The shifts are the result of 
a subtle balance of the repulsive Pauli blocking effect and attraction implied by the 
coupling of the D + meson to resonance- hole modes. 

§4. Conclusions 

The pion-mass dependence of the baryon octet and decuplet masses was dis- 
cussed. A novel approach was reviewed in which the latter are a solution of a set 
of coupled and non-linear algebraic equations. This is a direct consequence of self 
consistency imposed on the partial summation, i.e. the masses used in the loop 
functions are identical to those obtained from the baryon self energies. As a strik- 
ing consequence a discontinuous dependence of the baryon masses on the pion mass 
arises. Typically the baryon masses jump at pion masses as low as 300 MeV. Most 
spectacular is the behavior of the S mass. At small pion masses it decreases with 
increasing pion masses. At a critical pion mass of about 300-400 MeV it jumps up 
to a value amazingly close to the prediction of the MILC collaboration. 

In the second part of this talk we discussed recent results on the properties 
of D± and Df mesons in cold nuclear matter based on self consistent coupled- 
channel dynamics. It was pointed out that the D + ,D^ spectral distributions are 
strongly distorted in a nuclear medium due to the presence of resonances that couple 
strongly to the D + ,Djr nucleon channel. Since even the existence of most resonances 
that are predicted by coupled-channel dynamics is not yet confirmed experimentally 
such studies suffer from large uncertainties. This asks for a dedicated experimental 
programm to unravel the spectrum of baryons with charm content. 
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